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ABSTRACT 

We present Herschel observations at 70, 160, 250, 350 and 500 jiva of the environment 
of the radio galaxy 4C+41.17 at z = 3.792. About 65% of the extracted sources are 
securely identified with mid-IR sources observed with the Spitzer Space Telescope at 
3.6, 4.5, 5.8, 8 and 24 /im. We derive simple photometric redshifts, also including 
existing 850 /xm and 1200 fim data, using templates of AGN, starburst-dominated 
systems and evolved stellar populations. We find that most of the Herschel sources 
are foreground to the radio galaxy and therefore do not belong to a structure associated 
with 4C+41.17. We do, however, find that the SED of the closest (~ 25" offset) source 
to the radio galaxy is fully consistent with being at the same redshift as 4C+41.17. 
We show that finding such a bright source that close to the radio galaxy at the same 
redshift is a very unlikely event, making the environment of 4C+41.17 a special case. 
We demonstrate that multi- wavelength data, in particular on the Rayleigh- Jeans side 
of the spectral energy distribution, allow us to confirm or rule out the presence of 
protocluster candidates that were previously selected by single wavelength data sets. 

Key words: galaxies: individual: 4C+41.17 - galaxies: clusters: general - galaxies: 
high-redshift - techniques: photometric. 
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f Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 



1 INTRODUCTION 

1.1 High- Redshift Radio Galaxies as Tracers of 
Protoclusters 

High-redshift radio galaxies (HzRGs) are galaxies in the dis- 
tant universe [z > 1) showing enormous radio luminosities 



2 D. Wylezalek et al 



Stevens et all 120031: [Falder et all l20ld: 
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(Lsoomhz > KFWHz" 1 . iMilev & De Breuckl l2QQ8h . They 
are extremely rare objects, with numb er densities ^ 10~ 8 
Mpc~ 3 in the redshift range 2 < z < 5 jPunlop &; Peacockl 
Il99d : IWillott et aHl200ll : IVenemans et al.ll200/t ). Investigat- 
ing their spectral energy distribution (SED) reveals fea- 
tures of their stellar, dust and AGN components. In partic- 
ular, studies of the stellar and dust component have shown 
that HzRGs are a mongst the most massive galaxies in the 
early universe (e.g., Seymour et a l. 2007; B ryant et al ] l2009l : 
iDe Breuck et alJl201Ch/ According to the hierarchical model 
of galaxy assembly (| White & Reed [19781 ) , this implies that 
they reside in peaks of dark matter overdensities. As galaxy 
clusters represent the most massive structures in the uni- 
verse, HzRGs are expected to preferentially reside in sites of 
galaxy cluster formation. At z = 2 the universe is only ^3.2 
Gyr old and galaxy clusters are likely still forming but have 
not had time to virialize. For this reason we refer to these 
matter overdensities as protoclusters. Observations have in- 
deed shown that H zRGs preferentially reside in overdense 
environments (e.g. 
[ Stevens et alll2010l : 
I2012T ) and protoclusters are very likely to be found in the 
vicinity of these objects. As HzRGs are found up to very 
high redshift, they serve as efficient beacons for identifying 
very high redshift galaxy clusters. The fields of HzRGs are 
therefore unique laboratories to study the formation and 
evolution of the first galaxies and galaxy structures. 



1.2 The HeRGE Project 

With the launch of the Herschel satellite ([Pilbratt et al.l 
l2010h , it is possible for the first time to obtain full coverage 
of the far-IR SED for a large sample of HzRGs. The Her- 
schel Radio Galaxy Evolution project (HeRGE) makes use 
of the two imaging instruments on b oard Herschel: PACS , 
the Photodetecting Array Camera (|Poditsch et al] l20ld ) 
and SPIRE, the Sp ectral and Photometric Imaging Receiver 
(iGriffin et alJl20ld ). These instruments cover a wavelength 
range of 70 — 500/im and thus constrain the far-IR dust peak 
very well for a range of redshift s. The project was granted 
~ 27h of OT1 observing time (PI: N. Seymour) allowing 71 
HzRGs to be observed in five bands in PACS and SPIRE 
(PACS: 70/100 /xm, 160 /mi; SPIRE: 250 /mi, 350 /mi, 500 
/mi) . In additio n to studying the radio galaxies themselves 
in more detail |lvison et al.ll20l3 : iRocca-Volmerange et all 
l2012l : ISevmour et alJ2012h project HeRGE allows us, for the 
first time, to systematically study the environments of the 
radio galaxies at these wavelengths, reaching out 1-3' from 
the HzRGs. This complements our st atistical studies of the 
HzRG environme nts in the mid-IR ([Galametz et al.ll2012l : 
iMavo et alJl2Q12h . Reaching out to longer wavelengths al- 
lows us to constrain the dust peak of the SEDs and derive 
photometric redshift estimates to confirm or rule out over- 
densities associated with the HzRG. 

This work reports our pilot study of the well known HzRG 
4C+41.17. This analysis will be expanded systematically to 
the whole data set in future work. 

1.3 4C+41.17 

4C+41.17 at z = 3.792 is one of the best studied HzRGs. I t 
was discovered by Ch ambers, Milev, &; van Breugell (|l990h . 



The steep radio spectrum (a ~ —1.3) together with ex- 
tended optical continuum emission and the large rest frame 
Lya equivalent width (~ 270 A) identified 4C+41.17 as 



a HzRG. Its high far-infrared luminosity, Lf 



- mirared luminosi t y, -Lf ik ~ 1U i^q 
(iBenford et al.l Il999l : iHumphrev etafl 1201 lh . large dust 
mass (iDunlop et"alf Il994h and molecular gas reservoir 
(IDe Breuck et al.ll2005l " make this radio galaxy a very likely 
site of an enormous starburst at high redshift. Similar far- 
infrared luminosities h ave also been found for other high 
redsh ift radio galaxies ([Barthel et al.l 120121 : I Seymour et al.l 
l2012h accumulating the evidence for massive starburst in 
these galaxies. Deep observations at 450 /xm and 850 
/im carried out with SCUBA (Submillim ete r Common- 
User B olometer Array, I Holland et al.|[l999h bv llvison et al.l 
(l2000h in the field centered on 4C+41.17 show an order-of- 
magnitude overdensity of luminous sub-mm galaxies within 
a 2.5 arcmin diameter region centered on the radio galaxy. 
From tentative redshift constraints based on the 450 to 
850 /xm and the 850 /m i to 1.4 GHz flux density ratios 
of sources then available, llvison et al.l (|2000h conclude that 
the overdensity is consistent with lying at the same red- 
shift as the radio source, 4C+41.17, and therefore suggests 
a likely protocluster. However, photom etric redshifts esti - 
mated from the 1.6 fim stellar bump bv lGreve et all (|2007h 
place at leas t two o ut of the five sub-mm sources reported by 
llvison et al.l {[2000) at redshifts lower t han 1.3. Greve et al.l 
(l2007h also present deep SHARC-II ([Powell et al] l2Q03h 
350/xm and MAMBO (Ma x-Planck Millimeter Bolometer 
Array, iKrevsa et al.l [l998) 1200 /xm imaging of the field 
around 4C+41.17 and combine them with multi- wavelength 
data at 3.6, 4. 5, 5.8, 8 /mi from Spitzer IRAC (Infrared 
Array Camera, iFazio et al.ll200j ). 24 and 70 / mi data from 
Spitz er MIPS (Multiband Imaging Photometer, iRieke et al.l 
l2004h and 850 /xm observations from SCUBA. They find a 
surface density of ~ 0.24 1200 /mi sources per arcmin -2 to 
a depth of ~ 2 mJy, consistent with the blank field source 
density at this wavelength. From cross- correl ation analy- 
sis an d estimation of photometric redshifts, I Greve et al.l 
(l2007h conclude that at least half of the sub-mm galaxies 
are foreground sources and are not, in fact, associated with 
4C+41.17. 

In this paper, we present a multi- wavelength study of 
the environment of the HzRG 4C+41.17, recently ob- 
served within the HeRGE projec t in fiv e PACS and SPIRE 
bands. IRocca-Volmerange et all (|2012h present a full, in- 
depth study of the SED of the radio galaxy itself. We make 
use of the data at hand to derive photometric redshifts and 
to confirm or rule out the companionship of the galaxies in 
the field with the HzRG. Section 2 describes the observa- 
tions and reduction of the multi wavelength data. Section 3 
gives details of the source extraction and cross-correlation. 
In section 4 we present our analysis and draw conclusions 
in §5. Throughout the paper we assume Ho = 70 km s _1 

M P C-\ Onatter = 0-3, 0\ = 0.7. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Far-Infrared Observations 

Observations at 70 and 160 /xm were obtained with the Her- 
schel/PACS instrument on UT 2010 October 12. The image 



covers ~ 20 arcmin 2 . We retrieved the Level data from 
the Herschel Science Archive and processed it using version 
7.3.0 o f the Her schel Interactive Processing Environment, 
HIPE (|Ottll2010h . The data were taken to Level 1 following 
the standard pipelines provided in HIPE. To create Level 2 
products, we slightly adapted the standard pipeline to cor- 
rect for the slew to target data and to improve the point 
source sensitivity by decreasing the high pass filter radiufl 
The Herschel /SPIRE instrument observed a region cover- 
ing - 80 arcmin 2 around 4C+41.17 on UT 2010 Septem- 
ber 21 with all three bands, at 250 /im, 350 /im and 500 
/im. The exposure times for the PACS/SPIRE observations 
were 2 x 1404s and 721s, respectively, reaching an average 
la depth of 6.0, 6.4, 10.2, 9.6, 11.2 mJy at 70, 160, 250, 
350 and 500 /im, respectively. Both the SPIRE and PACS 
observations are part of the guaranteed time key program 
The Dusty Young Universe: Photometry and Spectroscopy 
of Quasars at z > 2 (Observation ID: 1342206336/7 and 
1342204958, PI: Meisenheimer) . 
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2.2 Mid-Infrared Data 

In addition, we include S pitzer IRAC and MIP S observa- 
tions in the analysis from ISevmour et al.l (|2QQ7h . The field 
was deeply mapped using all four IRAC bands (3.6, 4.5, 5.8, 
8 /im - referred to as channels 1, 2, 3 and 4), covering an 
area of 5.3 x 5.3 arcmin 2 , and all three MIPS bands (24, 70 
and 160 /im), covering an area of ~ 8.0 x 7.4 arcmin 2 . The 
exposure times were 5000s for the IRAC observations and 
267 s, 67 s and 2643 s for the three MIPS bands, in order 
of increasing wavelength. The data were reduced using the 
Spitzer reduction package, MOPEX. In this work, we only 
use the 24 /im images given the deeper PACS observations 
at longer wavelengths. The 3a depths reached were 0.8, 1.1, 
3.2 and 4.3 /iJy for the IRAC channels 1, 2, 3 a nd 4, respec- 
tively and 30 /iJy for the MIPS 24 /im image ([Greve et al.l 
l2007t ). 



2.3 (Sub)millimetre Data 

A field covering ~ 58 arcmin 2 around 4C+41.17 was imaged 
at 1200 /im with MAMBO. Details of t he observations, dat a 
reduction and analysis are reported by IGreve et al.l (|2007l ). 
Positions and flux densities of the extracted sources are 
taken from there. 

4C+41.17 was also observed at 850 /im with SCUBA (cov- 
ering an area ^ 2.5 arcmin in diameter, llvison et al. I l2000l : 



Stevens et al 



Iyisonetal 



.20031) . The data were initially published by 
and details can be found there. 



1 http://herschel.esac.esa.int/twiki/pub/Public/ 
PacsCalibrationWeb /bolopsfvl .01 .pdf 



Table 1. Flux densities and la uncertainties of sources with at least two Herschel detections 



Source 




£>4.5nm 


S5.8fxm 




*S'24/Ltm 


$70 fj,m 




$250/zm 








*Sl200/zm 










[A*Jy] 


[mjy] 


[mjy] 


[mjy] 


[mjy] 


[mjy] 


[mjy] 


[mjy] 


[mjy] 


1 
















54 ± 15 


61 ± 13 


42 ± 11 






2 


41 ±4 


46 ±5 


44 ±5 


51 ±5 


0.37 ±0.03 






47 ± 11 


38 ± 11 


36 ± 16 




3.0 ±0.6 


4 










0.48 ±0.03 






42 ± 13 


22 ±7 








5 


16 ±2 


20 ±2 


26 ±3 


13 ±2 


0.31 ±0.02 






19 ±3 


26 ±8 


10 ±3 




7.5 ±0.6 


7 


130 ± 13 


104 ± 10 


72 ±8 


193 ± 20 


0.39 ±0.03 






31 ± 13 


24 ± 11 


19 ±9 






9 


99 ± 10 


73 ±7 


79 ±9 


97 ± 10 


0.98 ±0.05 


19 ±5 


67 ±6 


64 ± 11 


31 ±6 








11 


29 ±3 


29 ±3 


24 ±3 


26 ±3 


0.35 ±0.03 




24 ±7 


44 ± 10 


39 ±8 


33 ±9 


9±1 


4.6 ±0.4 


12 














19 ±7 


27 ±9 


17 ±4 






3.0 ±0.6 


13 




27 ±3 




45 ±5 


0.32 ±0.02 






15 ±6 


30 ± 14 








16 


11 ± 1 


14 ±2 


24 ±3 


45 ±5 


0.47 ±0.03 




15 ±6 


42 ±5 


48 ±4 


39 ±4 


12 ±1 


3.8 ±0.4 


4C+41.17 


17 ±2 


20 ±2 


27 ±3 


31 ±4 


0.36 ±0.03 




16 ±7 


36 ±4 


43 ±4 


38 ±5 


12 ±1 


4.4 ±0.4 


18 














27 ±6 


33 ±8 


21 ±6 








19 


128 ± 13 


106 ±11 


98 ±11 


131 ± 13 


0.62 ±0.03 


8±3 


23 ±7 


29 ± 11 










21 


21 ±2 


26 ±3 


28 ±3 


22 ±2 








15 ±3 


36 ± 15 


21 ±9 




2.6 ±0.6 


24 


38 ±4 




51 ±5 










10 ±4 


11 ±5 






3.6 ±0.6 


28 










0.35 ±0.03 






37 ±4 


22 ±8 








29 


130 ± 13 


137 ± 14 


117 ± 12 


368 ± 37 


1.41 ±0.07 






13 ±4 


17 ±5 
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3 SOURCE EXTRACTION AND 

CROSS-CORRELATION ANALYSIS 

3.1 Source Extraction 

3.1.1 PA CS/ SPIRE Source Extraction 

Source extraction in the PACS and SPIRE images is per- 
formed using the tool sourceExtractorDaophot that is in- 
cluded in HIPE. The FWHM of the different bands are taken 
from the PACS Observer's Manual and are 5.2, 12, 18.1, 
25.2 and 36.3" for 70, 160, 250, 350 and 500 /im, respectively. 
The parameters for the source extraction, such as shape pa- 
rameters roundness and sharpness, are tuned such that false 
detection rates and source blending is minimized (Table [2]). 
We extract sources at a significance ^ 2.5cr within a circle 
of 3.3' (corresponding to 34.2 arcmin 2 ) radius around the 
radio galaxy. Due to the scanning mode the coverage is in- 
homogenous further away from the image center. We extract 
two sources from the PACS 70 /im image, 8 sources from 
the PACS 160 /im image, 27 sources from the SPIRE 250 
/im image, 16 sources from the SPIRE 350 /tm image and 8 
sources from the SPIRE 500 /im image. The extracted source 
positions and given names are listed in Table [3] in order of 
increasing RA. We derive aperture photometry for the ex- 
tracted sources applying an aperture correction of 1.45 and 
1.44 to the blue and red PACS flux densities, respect ivefjfl 
Due to the inhomogeneous coverage in the Herschel images 
the uncertainty on the flux densities is derived from sky an- 
nuli (see Table [2]) around each source. 

Aperture photometry is, however, not applicable in the case 
of source 16 and 4C+41.17, which are blended. We therefore 
apply PSF photome try to those sources using StarFinder 
(jDiolaiti et alJl20QCh . a code designed to analyze images in 
very crowded fields. The deblending strategy in StarFinder 
consists of an iterative search for residuals around the object 
and subsequent fitting. We assumed the PSF to be Gaus- 
sian with a FWHM corresponding to the beam size. The 
positions of the two sources were determined independently 
in each Herschel image as different material is probed at 
different wavelengths. The flux density measurements with 
StarFinder are consistent with the ones obtained with HIPE 
for unblended sources. Postage stamps of the image, syn- 
thetic image and residual image after deblending is shown 
in Figure [TJ No other sources in the field are blended and 
the fluxes densities and uncertainties (including the 15% and 
7% flux calibration uncertainties added in quadrature to the 
statistical un certainties for PACS and SPIRE flux densities, 
respectively, ISevmour et~aT1l2012h are given in Table [1] 



3.1.2 IRA C/MIPS Source Extraction 

Source extraction is performed using SExtractor 
(Ber tin fc Arnoutd Il996h in dual image mode using 
the 4.5 /xm image for detection. We only re port sources 
detect ed with a significance ^ 3a. Unlike iGreve et al.l 
(l2007h . we use a smaller, 4" diameter aperture for the 
IRAC images because of the close proximity of other 
sources in that crowded field. Tests with an aperture of 



2 http://www.iac.es/proyecto/herschel/pacs/pacs_om.pdf, p. 13 

3 http:/ /herschel. esac.esa.int/Docs/PACS/html/pacs_om. html 
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Figure 1. 0.6 ; x 0.6' postage stamps of the data (left), synthetic 
image derived by StarFinder (center) and residuals (right) for 
source 16 and 4C+41.17. From top to bottom the images at 250, 
350 and 500 /im are shown. The sources are blended in all three 
images but the homogeneous residual image shows the good de- 
blending with StarFinder. Red crosses indicate the positions of 
4C+41.17 (upper source) and source 16 (lower source). 



9.26" diameter show that flux from neighbouring sources 
results in overestimated flux densities (e.g. for 4C+41.17 
itself). We apply aperture corrections of 1.205, 1.221, 1.363, 
1.571 to IRAC channels 1, 2, 3 and 4, respectively. MIPS 24 
/mi flux densities are measured in 5.25" aperture radii. The 
aperture correction applied, 1.78, is calculated as described 
by the MIPS instrument handbookQ. The uncertainties 
reported in Table [1] include the 10% and 4.5% systematic 
uncertainties for the IRAC and MIPS flux densities, re- 
spectively, that were added in quadrature to the statistical 
uncertainties to account for the absolute flux calibra tion 
and color correction uncertainties ([Seymour et al.ll2007h . 

3.2 Cross-Correlation Between Bands 

After extracting sources in the different images with very 
different spatial resolutions we cross- correlate the sources 
in order to derive a clean, mult i- wavelength source catalog. 
We only consider the 17 sources that have at least two de- 
tections in the Herschel bands in order to minimize false 
detections. For the cross-correlation, we choose the SPIRE 
250 /im whose la positional accuracy (~ 0.6 x jf[ out- 

performs the other bands. Although the PACS images have 
an even better spatial resolution, they cannot be used sys- 
tematically as reference images due to their shallowness and 
small field of view (see Figure [2]). We then cross- correlate 
the cleaned source list with the sources detected at shorter 

4 http: / / irsa. ipac.caltech.edu / data/SPITZER/docs / mips / 
mipsinstrument handbook /50/ 

5 http:/ /herschel. esac.esa.int/hcss-doc- 

8. 0/load/hcss_urm/html/herschel.ia. toolbox, srcext. 
SourceExtractorDaophot Task, html 
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Table 2. HIPE parameters used to extract sources with sourceExtractorDaophot and details of the observations in different Herschel 
bands. Roundness and sharpness parameters between -2 to 2 and -1.5 to 2 were used for the extraction in all bands. 



Band 


FWHM 


beam area 


aperture radius 


sky annulus 


exposure time 


average la depth 




[arcsec] 


[arcsec 2 ] 


[arcsec] 


[arcsec] 




[mJy] 


PACS 70 yum 


5.2 


30.6 


6 


6-10 


2 x 1404s 


6.0 


PACS 160 lira 


12 


163.2 


10 


10-15 


2 x 1404s 


6.4 


SPIRE 250 lira 


18.1 


373.3 


22 


22-32 


721s 


10.2 


SPIRE 350 /im 


25.2 


716.7 


30 


30-40 


721s 


9.6 


SPIRE 500 /im 


36.3 


1493.1 


42 


42-52 


721s 


11.2 



Table 3. Astrometry of the Herschel sources in the 4C+41.17 field, listed in order of increasing RA. Positions are in the J2000 system. 



Source 


RA 70^,m 


Dec 


•70^,m 


RA 160/nm 


Dec -160Aim 


RA 250Aim 


Dec 


250/xm 


RA 350Atm 


Dec 


350Aim 


RA 500Atm 


Dec 


500Atm 


1 












6:50:37.0 


41 


28:54 


6:50:36.9 


41 


28:53 


6:50:36.6 


41 


28:44 


2 












6:50:40.5 


41 


30:03 


6:50:40.7 


41 


30:05 


6:50:40.7 


41 


30:05 














6:50:40.9 


41 


32:35 














4 












6:50:42.6 


41 


28:27 


6:50:42.8 


41 


28:25 








5 












6:50:43.1 
6:50:44.0 


41 
41 


29:21 
27:43 


6:50:43.7 


41 


29:16 


6:50:43.2 


41 


29:17 


7 












6:50:45.6 
6:50:47.1 


41 
41 


32:40 
27:37 


6:50:44.9 


41 


32:37 


6:50:46.4 


41 


32:52 


9 


6:50:47.4 


41 


30:46 


6:50:51.1 


41:30:06 


6:50:47.2 


41 


30:44 


6:50:47.3 


41 


30:46 








10 












6:50:47.4 


41 


28:11 














11 








6:50:48.9 


41:31:27 


6:50:48.8 


41 


31:26 


6:50:49.1 


41 


31:28 


6:50:49.0 


41 


31:22 


12 








6:50:50.1 


41:28:21 


6:50 


50.1 


41 


28:19 


6:50:50.0 


41 


28:21 








13 












6:50 


50.3 


41 


33:01 


6:50:50.5 


41 


32:58 








14 












6:50 


50.3 


41 


27:37 














15 












6:50 


51.2 


41 


29:48 














16 








6:50:51.2 


41:30:06 


6:50 


51.4 


41 


30:07 


6:50:50.9 


41 


30:04 


6:50:51.1 


41 


30:06 


4C + 41.17 








6:50:51.9 


41:30:32 


6:50 


51.7 


41 


30:32 


6:50:52.1 


41 


30:31 








18 








6:50:52.4 


41:28:53 


6:50 


52.2 


41 


28:51 


6:50:52.4 


41 


28:55 








19 


6:50:54.4 


41 


29:33 


6:50:54.4 


41:29:33 


6:50 


54.2 


41 


29:34 














20 








6:50:54.6 


41:30:47 






















21 












6:50 


54.8 


41 


32:36 


6:50:54.8 


41 


32:29 


6:50:53.8 


41 


32:42 


22 


























6:50:56.1 


41 


33:35 


23 












6:50 


56.4 


41 
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and longer wavelengths. 

We look for MIPS counterparts within 10" of the 250 /im 
sources which corresponds to about the half-width at half 
maximum for the 250 /im observations and which also cor- 
responds to their 3a positional error assuming the bulk of 
th e 250 /im detections h as a signal to noise ratio (SN) of ^ 
3 dMagnelli et al.l l2Q12h . Following Sutherland & Saundersl 
(|l992h , we calculate the reliability R — exp(—7rr 2 <7i<72N) 
of finding no random source closer than the nearest can- 
didate where N is the number density of background ob- 
jects, r = \J (d\/ai) 2 + (d^/cr?) 2 is the normalized distance, 
di and cfo are the positional differences in each axis be- 
tween the sources, and o\ and 02 are the standard devi- 
ations of the error ellipse. Since <7\ — a 2 for our case, R 
simplifies to R = exp{-^{d\ + dl)N). iMavo et al.l (|2012h 
find a surface density 0.549 arcmin -2 at the depth of the 
M IPS image. This is in agreement with the density found 
bv lPapovich et all (|2004h for the 23.1 /iJy depth of the MIPS 
image of 4C+41.17. We adopt this value for calculating the 
reliability of MIPS-250 /im sources counterparts. Candidates 
with reliabilities R ^ 90% are counted as correct identi- 
fications. We then cross-correlate the MIPS identifications 
with the IRAC catalogs and compute the reliability in the 
same way using a su rface density of 2.80 arcmin -2 from 



iGalametz et all (|2012). For sources with no MIPS identifi- 
cation, we cross-correlated the 250 /im sources directly with 
the IRAC catalogs. Again, candidates with reliabilities ^ 
90% are counted as correct identifications. Although we re- 
quire R ^ 90%, the reliabilities for MIPS identifications are 



all above 96% and for IRAC identifications above 98%. Ta- 
ble 31 lists the separations between the MIPS and 250 /im 
counterparts and the IRAC and MIPS counterparts together 
with the corresponding reliabilities of them to be the cor- 
rect identifications. Following this methodology, out of the 
17 Herschel 250 /im sources, we identified MIPS counter- 
parts for 11 sources (65%), and IRAC counterparts for an 
overlapping, but not identical set of 11 far-IR sources (65%). 
We checked the identifications where better resolved PACS, 
intermediate wavelength data were available and confirmed 
our identified sources. 

In the Appendix we show postage stamps of all sources 
with SPIRE 250 /im sources that have high probability (R 
> 90%) counterparts in at least two IRAC bands and are 
detected in at least two Herschel bands. We also overplot 
the sub-mm position from iGreve et al ] (120071) where avail- 
able. Except for sources 19 and 24, the IGreve et al.l (l200<t ) 
astrometry is in very good agreement with our mid-IR iden- 
tifications. 



4 ANALYSIS 

4.1 Photometric Redshifts 

We now derive photometric redshifts in order to investigate 
a physical connection between the radio galaxy and the ob- 
jects in its vicinity. The following issues must be kept in 
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Table 4. Separation between MIPS and 250 /im sources, between 
IRAC and MIPS sources, and calculated reliabilities for the near- 
est source to be the correct identification. Where no value is given, 
the IRAC and/or MIPS source is either outside the field of view 
or there are no detections within a 10" radius. Italized numbers 
indicate that there are only IRAC2 and 4 or IRAC1 and 3 detec- 
tions for the 250 jum source as the other bands are outside the 
field of view. In case of source 21 the IRAC-250 //m position and 
corresponding probability is listed. 



Source 


MIPS-250/im 
separation 
["] 




IRAC-MIPS 
separation 
["] 


^IRAC 


1 

2 


9.18 


0.96 


2.81 


0.98 


5 


4.37 


0.99 


0.58 


0.99 


7 


3.70 


0.99 


0.89 


0.99 


9 


4.99 


0.99 


0.62 


0.99 


11 


5.01 


0.99 


2.46 


0.99 


12 










13 


2.54 


0.99 


0.51 


0.99 


16 


0.60 


0.99 


0.59 


0.99 


4C+41.17 


6.46 


0.98 


0.96 


0.99 


19 


4.86 


0.99 


0.71 


0.99 


21 






2.47 


0.99 


28 


4.87 


0.99 






29 


5.26 


0.99 


0.40 


0.99 



mind when deriving photometric redshifts from combined 
sub-mm, far-IR and near-IR observations: 

(1) As the far-IR emission is of thermal origin, chang- 
ing the dust temperature has the same effect on the 
sub-mm /mm colo rs as shifting the spectrum in redshift 
(Bla in et al.ll2002l ). It is thus impossible to estimate the red- 
shift from far-IR data alone; supporting observations are 
necessary to constrain a redshift. 

(2) The dust (at far-IR, sub-mm wavelengths) to stel- 
lar flux density (at near-IR wavelengths) ratio has a range 
of about 3 decades and varies with morphology, total IR 
luminosity and gas-phase metallicity ([Skibba et al.l 1201 lh . 
This is not always well represented in the available tem- 
plate libraries. When deriving our own empirical templates, 
we therefore create templates with a wide range of dust to 
stellar ratios, ranging from 100 to 5000. 

(3) Differing error bars for near-IR and far-IR observa- 
tions will introduce a bias in the photometric redshift fitting 
procedure, giving the high signal-to-noise ratio of the IRAC 
data more weight. 

For all 11 sources with detections in more than five wave- 
length bands we calculate photometric redshif ts using the 
code hyperz (iBolzonella. Miralles. & Pellcl2000h which min- 
imizes the reduced x? to find the best photometric red- 
shift solution. We use both synthetic and empirical AGN 
and starburst templ ates from the SWIRE template library 
(jPolletta et al.ll2007h complemented with our own newly de- 
rived templates. The latter are obtained by combining a 1 
Gyr old stellar popul ation template from the Pegase.2 spec- 
tral evolution model (|Fioc &; Rocca-Volmerangell 19991 ) dom- 
inating the near-IR emission and empirical dust templates 
dominating the far-IR/sub-mm emission. Three dust tem- 
plates are derived by (1) fitting the dust peak of 4C+41.17, a 



typical AGN dominated galaxy at high redshift; (2) the dust 
peak o f the lensed "eyelash" galaxy at z = 2.3 (S MM J2135- 
0102, llvison et alll20ld : ISwinbank et al.l l2Q10h , a typical 
starburst galaxy at high redshift; and (3) source 11, which 
is very well sampled at far-IR/sub-mm wavelengths and for 
which the spectr oscopic redshift is known (z sp ec = 1.18, 
iGreve et al.l 120071 ). In this way, not only templates derived 
from lower redshift galaxies, such as the SWIRE templates, 
are available to us but also templates derived from higher 
redshift galaxies. For each of the three dust templates we 
create 15 composite templates with different ratios between 
the stellar emission in the near-IR and the dust emission 
in the far-IR. In order to get a matching wavelength cover- 
age of the self-derived templates and the SWIRE templates 
we extend our templates by using greybody fitting results 
(see Section 4.2) for A > 1200 /xm. We then extract the 
best-fitting templates from our 45 self-derived templates and 
the SWIRE library templates. Ultimately, a set of 9 differ- 
ent templates (see Table [5} are used for the 11 sources for 
which we derive photometric redshifts. Templates 1 (Spiral 
C) and 3 (starburst) (both from the SWIRE library) are 
gener ated from the SEP of these objects using the GRASIL 
code (jSilva et al.lll998T ) and improved by using IR spectra 
from the PHT-S spectrometer on the Infrared Space Obser- 
vatory and from IRS on Spitzer (|Houck et al.l 12004 ). Tem- 
plate 2 is an empirical composite AGN+starburst template 
that fits IRAS 19254-7245. Template 4, 5, 6, 7, 8 and 9 are 
new, self-derived templates, with their properties described 
in Table [5] The resulting x 2 distribution and the best x 2 are 
thus derived by considering all redshifts and all templates 
in the final set. Note that the final x 2 curve shows the min- 
imum x 2 f° r the template set as a function of redshift and 
therefore is dependent on the template set used. 
Because of varying spatial coverage of the multi- wavelength 
data, filters are ignored for 'out-of- field' sources, but when 
a source is observed, but undetected, 3a upper limits are 
taken into account by hyperz. We present the results of our 
photometric redshift estimates in Table [6] and show the best- 
fitting SEDs in the Appendix. 

As mentioned above, differing errors bars for the near-IR 
and far-IR observations introduce a bias in the fitting proce- 
dure giving the high signal-to- noise IRAC data more weight. 
However, by allowing a range of various ratios between the 
stellar (near-IR) and dust (far-IR) emission in the fitting 
templates we already make sure that the fits are not depen- 
dent on the IRAC data only but that also the relative con- 
tribution of the sources of emission is taken into account. To 
test this in more detail, we repeat the fitting procedure by re- 
laxing the IRAC uncertainties to 20%. The best-fit redshifts 
are in agreement with the previously derived ones within the 
uncertainties. This shows that our results are not strongly 
biased by th e IRA C data. 

IGreve et al.l (l2007h derived photometric redshifts using the 
1.6 fim rest-frame stellar 'bump' in the observed IRAC data. 
They also estimate the redshift from the radio/sub-mm/mm 
color, but these only yield crude estimates, consistent with 
the redshift estimation from the stellar bump. Therefore, we 
only list Zbump in Table [6] which compares to our photomet- 
ric redshifts and spectroscopic redshifts that exist for some 
of the sources. The uncertainties listed in Table[6]only reflect 
the la formal uncertainties near the minimum of the x 2 dis- 
tribution and may be severely underestimated. We describe 
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Table 5. Summary of the templates used for deriving photometric redshifts. Template 1-3 are from the SWIRE template library 
(IPolletta et al.l [2007 N ) . Template 4-9 are newly derived by combining the far-IR emission of 4C+41.17, the eyelash galaxy (SMM J2135- 
0102) and source 11, for which the spectroscopic redshift is known, with a 1 Gyr old stellar population from the Pegase.2 spectral 
evolution model (|Fioc fc Rocca-Volmerangel 11999). The far-IR and stellar emission were normalized to their peak flux densities and 
combined with varying ratios, as indicated. 



Template ID Description 



1 Spiral C galaxy template, SWIRE template library 

2 Seyfert 2+Starburst/ULIRG template for IRAS 19254-7245, SWIRE template library 

3 Starburst/ULIRG template for IRAS 20551-4250, SWIRE template library 

4 4C+41.17 far-IR template + old stellar population, stellar peak to dust peak ratio: 300 : 1 

5 4C+41.17 far-IR template + old stellar population, stellar peak to dust peak ratio: 700 : 1 

6 4C+41.17 far-IR template + old stellar population, stellar peak to dust peak ratio: 4500 : 1 

7 SMM J2135-0102 far-IR template + old stellar population, stellar peak to dust peak ratio: 500 : 1 

8 source 11 far-IR template + old stellar population, stellar peak to dust peakratio: 700 : 1 

9 source 11 far-IR template + old stellar population, stellar peak to dust peak ratio: 1700 : 1 



Table 6. Photometric redshifts, Zphot? were derived with hyperz 
for all sources with at least five detections. Stellar bump pho- 
tometric redshifts, ^j^mp, and spectroscopic redshifts, 

2-spec are 

from I Gr eve et all ([20071 ). Template ID's are described in Table [5] 



Source 


^phot 


template ID 


^bump 


^spec 


2 


2.5±0.4 


3 


~ 1.8 




5 


2.4±0.2 


3 


~ 2.6 


2.672± 0.001 


7 


0.5±0.1 


1 






9 


0.6±0.1 


8 






11 


1.2±0.2 


9 


< 1.3 


1.184 ± 0.002 


13 


2.2±0.4 


5 






16 


4.0±0.1 


6 






4C+41.17 


3.5±0.2 


2 


- 4 


3.792 ± 0.001 


19 


2.0±0.1 


4 


< 1.3 


0.507 ± 0.020 


21 


2.7±0.2 


7 


~ 1.8 




29 


l.OdtO.l 


4 







the quality of the photometric redshifts for each source in 
the Appendix. Figure [2] shows the spatial distribution of the 
sources, with white, open stars representing lower redshift 
objects (z < 3) and filled stars representing objects with 
z > 3. 

The average redshift for these 11 sources is 2.0 =b 0.8. 
This agrees with the average redshift for SPIRE 250 
jim. selected sources, z — 1.8 =b 0.2, recently found by 
iMitchell- Wynne et all (120121 ). 

Most of the SPIRE-selected sources are found to be at 
z < 2.5, ruling out any physical connection with the ra- 
dio galaxy and confirming that most of the far-IR sources 
in the vicinity of 4C+41.17 are likely foreground. Only one 
source, object 16, potentially lies at the same redshift as 
4C+41.7. The x 2 distribution of this source shows a clear 
dip at z = 3.8. We therefore assume that object 16 and 
4C+41.17 are at the same redshift, z = 3.8, and adopt this 
assumption for our subsequent analysis. 



4.2 Herschel non-detections 



llvison et all (|2010h report five 850 fim detected sources, of 
which two are described as marginal detections; for the ro- 
bust subset, we find viable Herschel counterparts to two. We 
find the sa me detection rate for the robust 1200-/im sources 
reported by iGreve et al.l (|2007h . The fact that they are not 



o 

LU 
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A R.A. [arcmin] 

Figure 2. Coverage map and spatial distribution of sources with 
derived photometric redshifts, centered on 4C+41.17. Dark pix- 
els indicate regions with low coverage. White, open stars indicate 
sources that have z p hot < 3, orange, filled stars show sources 
with Zphot > 3. 4C+41.17 and source 16 are very likely to be at 
the same redshift. The dashed circle (r = S.S^/ellipse (a = 2.5', 
b = 1 ; ) indicates the area for identifying targets in the SPIRE and 
PACS maps, respectively. The dotted box shows the coverage of 
IRAC1 and IRAC3, the dashed-dotted box shows the coverage of 
IRAC2 and IRAC4 and the long dashed circle shows the cover- 
age of SCUBA. The MIPS and MAMBO images cover the whole 
SPIRE area and are therefore not illustrated here. The cover- 
age/error in the region used for extraction is not homogeneous. 



detected in any other wavelength band may suggest some 
of them are just statistical fluctuations. This would be es- 
pecially true for those that are only marginally detected. 
The 1200 /im sources, however, are all observed at a sig- 
nificance ^ 5cr. If the sources are real they are likely very 
dust obscured sources belonging to the high redshift (z > 4) 
tail of sub-mm bright sta r- forming galaxies ([Swinbank et al I 
12012b IWalter et alll2012h . None of these sources has signifi- 
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Table 7. Derived dust temperatures (T^), grain emissivity indices 
(P), far-IR luminosities (Lfir) and star formation rates (SFRs) 
for sources with at least 4 detections in SPIRE, MAMBO and 
SCUBA. Spectroscopic redshifts where used where available and 
are marked in italics. For source 16, we assumed the redshift of 
4C+41.17 (z = 3.8). 



Source 


z 


T d 


P 


^FIR 


SFR 






[K] 




[10 13 Lo] 


[M^yr- 1 ] 


2 


2.6 


38 


1.6 


0.7 


1200 


5 


2.7 


40 


1.5 


0.3 


500 


11 


1.2 


26 


1.4 


0.1 


200 


16 


3.8 


48 


1.7 


1.8 


3100 


4C+41.17 


3.8 


52 


1.6 


1.6 


2800 


21 


2.7 


32 


1.7 


0.4 


700 



cant MIPS detections or unambiguous IRAC detections and 
we are therefore not able to estimate a likely redshift range. 



4.3 Far-IR Luminosities, Star-Formation Rates 
and Limits 

For sources with more than 4 detections in the the SPIRE, 
MAMBO and SCUBA bands we derive dust temperatures, 
far-IR luminosities and star formation rates (SFR). Sources 
2, 5, 11, 16, 4C+41.17 and 21 were fitted with a grey-body 
law of the form: S u oc 'B U (T) = - h "/ kTd — ^ where S u is 
the flux density at the rest-frame frequency u, f3 the grain 
emissivity index and Td the dust temperature. Dust tem- 
peratures for an interstellar medium only heated by star- 
formation in expecte d to range b etween ~ 20-60 K; ft can 
range between 1-2.5 (|Casevll2012h . Far-IR luminosities were 
derived by integrating their SED over the wavelength range 
40-1000 /im and applying the relation Lfir = AttD^Ffik 
where Dl is the luminosity distance computed from their 
photometric redshifts. Where spectroscopic redshifts were 
available those were applied. Source 16 was assumed to be 
at the redshift of the radio galaxy (z = 3.8). We then es- 
timated the star fo rmation rates b y using SFR [Mq] = 
Lfir/5.8 x 10 9 L (lKennicuttlll998h . The results are listed 
in Table □ 

Given the shallowness of the SPIRE images, at z = 3.8, 
we are only sensitive to the most massive starbursts. As- 
suming a dust emission from the starburst with P = 1.5 
and Td = 45 K, we find that at z = 3.8 we can only detect 
galaxies with a SFR > 2600 MQyr -1 . We therefore can only 
report on the presence of strongly starbursting galaxies in 
the field of 4C+41.17. 



4.4 Number Density 

In order compare the source d ensity to the Her schel wide 
field s urveys Herschel- A TLAS (jEales et al.ll2010h and Her- 
MES ([Oliver et al.ll2012h , we restrict this analysis to a flux 
density limit at which our catalogs are complete. We esti- 
mate the incompleteness by placing artificial sources in our 
images and applying the source extraction algorithm on the 
modified images. The number of successful recoveries then 
provides us with an estimate of the incompleteness for var- 
ious flux density bins. A completeness of ~ 95% is reached 
at ~ 35m Jy (corresponding to a SFR ~ 2600 MQyr -1 at 



z = 3.8) for our 250 /am catalog. We find 8 sources above 
that flux density limit in the extraction area of 34.2 arcmin 2 , 
resulting in a surface density of ~ 0.23 =b 0.08 arcmin -2 for 
the field of 4C+41.17. The SPIR E blank field number co unts 
at our flux density limit found by Cl ements et al 1 (|2010h are 
N(S25Qtxm > 36 mJy) - 0.121 ± 0.002 arcmin -2 suggest- 
ing a marginally signi ficant overde nsity of a factor ~ 2. On 
the other hand, [Oliver et al 1 (|2010h find blank field number 
counts of N(S25Qnm > 30 mJy) ~ 0.18 ± 0.15 arcmin -2 . 
We find 10 sources with flux densities above 30m Jy corre- 
sponding to a surface density of 0.3 ±0.1 arcmin -2 . As our 
catalogs are not yet complete at that limit this result has to 
be treated with caution but still hints to a slight overdensity 
for the field of 4C+ 41.17. 

iMavo et all (|2012h find a density of MIPS 24 fim sources 
which is ~ 2 times higher than the blank field mean density, 
which agrees with our Herschel observations. Nevertheless, 
this field was not counted as significantly overdense as it is 
still less than 3a above the mean. Compared t o the mean 
density of HzRG fields analyzed by IMavo et al.l ([20121 ) , the 
field around 4C+41.17 is typical as compared to other ra- 
dio galaxies in terms of density (£40+41.17 = 2.1 arcmin -2 , 
compared to (fazRG) ~ 2.2 arcmin -2 ). 

Selecting color-s elected high-redshift IR AC sources in the 
fields of HzRGs, Gal ametz et al ] (l2012h finds 4C+41.17 to 
be overdense with a compact clump of IRAC sources identi- 
fied r>j Y south of the radio galaxy. However, the IRAC color 
criterion applied simply identifies sources at z > 1.2. Given 
that we find an exces s of galaxies at z ~ 2.5 (see Figure [2} , 
the clump detected bv lGalametz et ah l (|2012h is likel y afore- 
ground structure, as also suggested by Greve et al.l (|2007h . 
Considering all data at hand, we therefore find no indica- 
tions for a remarkable overdensity in the field of 4C+41.17. 
In the 250 /xm image, source 16 is found ~ 25 arcsec from 
the 250 fim position of 4C+41.17 (~ 180 kpc at z = 3.8). 
The probability of finding such a bright far-IR source (42 
mJy at 250 /im, 48 mJy at 350 fim) at this distance to 
the H zRG is 4% and 8% for 250 /xm and 350 /xm, respec- 
tively ([Oliver et al 1l2010h . This probability is not remarkable 
when considering the whole sky but is very special when the 
evidence points to the two sources lying at the same red- 
shift. The probability of findi ng a 250 /im selected sourc e 
at z = 3.8 is also only - 5% ([Mitchell- Wynne et al.ll2012h . 
The probability of finding a source of that flux density at 
25" distance at z = 3.8 will therefore be <C 4% and is a very 
unlikely event. 



5 SUMMARY AND CONCLUSIONS 

Using Herschel PACS and SPIRE observations combined 
with Spitzer mid-IR observations, we have carried out a 
mult i- wavelength study of the environment of 4C+41.17, 
a powerful radio galaxy at z = 3.8. This pilot study for 
the HeRGE project clearly shows that far-IR observations 
combined with shorter wavelength observations improve our 
ability to securely disting uish overdensities found by differ- 
ent s election criteria (e.g. iGalametz et al.ll2012l : IMavo et al.l 
l2012h from truly clustered structures. The field of 4C+41.17 
has long been thought to host a galaxy c l uster associ- 
ated w ith the radio galaxy ([ivison et al 1 l200(t ). lGreve et al.l 
(2007|) already concluded from stellar bump photometric 
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redshifts that most of the sources might belo ng to a fore- 
groun d structure. Only source 16 (J065051.4 in iGreve et al] 
120071 ) appeared possibly associated with 4C+41.17. In this 
work we find strong indications that these two sources lie 
at the same redshift a nd thus that t here is a p hysical asso- 
ciation between them. Ilvison et al.l (|2000l ) and IGreve et al.l 
(2007) find that the radio galaxy and source 16 appear to 
be connected by a faint bridge in both the SCUBA and 
MAMBO map increasing the likelihood that this source is 
part of the same system as 4C+41.17. Source 16 makes the 
environment of 4C+41.17 special as the probability of find- 
ing such a bright source that close (~ 25 arcsec distant, 180 
kpc at z = 3.8) is only ~ 5%. 

However, close- by companion s ources might actua lly be 
a common feature for HzRGs. Ilvison et all (l2008h finds 
two clumps of emission 3.3 arcsec distant from the HzRG 
4C+60.07 tha t are most l i kely m erging with the z = 3.8 
radio galaxy. Ilvison et al.l (|2012h finds a bright sub- mm 
galaxy near the radio galaxy 6C 1909+72 that is most 
likely sha ring the same node o r filament of the cosmic 
web. Also. lNesvadba et al.l (|2009h find two CO-emission line 
components at a distance of ~ 80 kpc from the HzRG 
TXS0828+193 [z = 2.6) which may be associated with a 
gas-rich, low-mass satellite galaxy. Although these compan- 
ions are found much closer to the HzRG than source 16 
is to 4C+41.17, these observations suggest that companion 
sources around H zRGs may be a common feature (see also 
Ilvison et al.ll20l3 ) . We find that most of the Herschel far-IR 
sources in the vicinity of 4C+41.17 are foreground sources. 
However, this does not rule out the presence of a cluster 
around 4C+41.17 as our observations are only sensitive to 
galaxies with SFRs > 2600 MQyr" 1 . 

Caution is needed when identifying overdensities from a sin- 
gle wavelength data set. With IRAC and MIPS data avail- 
able for all sources being observed by the HeRGE project we 
will be able to identify likely protocluster candidates around 
the HzRGs. However, 850 /im data are required to con- 
strain the Ray leigh- Jeans part of the SED. We have there- 
fore started a systematic SCUBA- 2 follow up campaign to 
map the full SPIRE area of the HeRGE fields. 



APPENDIX A: NOTES ON INDIVIDUAL 
SOURCES 

In the Appendix we give more details on all sources with 
photometric redshifts derived in this work. For each source 
60" x 60" IRAC and MIPS grey-scale postage stamps are 
shown in the first row. The second row shows PACS and 
SPIRE grey scale postage stamps, 100" on a side. The 
postage stamps are centered on the 250 /im centroid, in- 
dicated by the orange square. The blue circle in the upper 
row shows the 10" search radius for the cross correlation 
analysis, the smaller, red circle in all stamps shows the aper- 
tures for the IRAC, MIPS, PACS and SPIRE images. The 
green diamond indicates the MAMBO position (in case of no 
MAM BO position, the SCUBA position) from IGreve et all 
(|2007h . We also show the SEDs and the minimum % 2 as a 
function of redshift from hyper z. Black data points are mea- 
sured values, green arrows upper limits and red dashed lines 
the best-fitting redshifted template. Detailed notes for each 
source are given below. 
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Figure A2. The IRAC photometry of this source shows a very prominent stellar bump that is well fit by the starburst template, leaving, 
no doubt on the low redshift ( £ n hot ~ 2.4) of this source. This is also confirmed by its spectroscopic redshift (z sp ec = 2.672) consistent 
with our photometric redshift. I Greve eta D (120071 ) finds a very extended Lycn halo extending 50 kpc from this source. 
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Figure A3. Although the far-IR photometric observations are not well fit by hyperz, the upper limit at 850/im is very constraining and 
places the dust peak at ~ 200 /im. The IRAC and MIPS photometric points can only be fit with a spiral template. The x 2 distribution 
places the source unambiguously at low redshift (z p hot — 0.5). 




Figure A4. The hyperz fit (source 11 dust template + old stellar population) fits the far-IR emission very well. The IRAC observations 
show the long wavelength tail of the stellar bump (at 1.6 /im restframe wavelength) indicating a low redshift. We find no secondary 
prominent dips in the % 2 distribution and this source is thus found to be at lower redshift (z p hot — 0-6) than the radio galaxy. 
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Figure A5. This source is detected in all far-IR bands (160-1200 /zm). The best % 2 is found with template 9 (this source + old stellar 
population) and gives a redshift of 1.2, matching the spectroscopic redshift of 1.184. 




Figure A6. Only five photometric points are available for fitting the SED of this source. The increasing emission towards 350 jim is 
not well fit and the best x 2 solution gives a redshift of 2.2. Blending of other sources in close proximity may cause this increasing flux 
at 350 i^m. 
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Figure A7. This source is nicely fit by an AGN dominated template, similar to the SED of 4C+41.17. The x 2 distribution shows a 
clear and prominent dip at the redshift of the radio galaxy. This source is therefore our most likely candidate to be associated with the 
radio galaxy. 
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Figure A8. The redshift of the radio galaxy is well constrained by the photometric redshift fitting using a composite AGN+starburst 
template (119254). A single significant dip appears at a redshift of ~ 3.5 in the % 2 distribution which is consistent with the spectroscopic 
redshift of 3.792. Note that this source is not fit by its own template as the stellar to dust peak ratio in those templates is not as optimal 
as in template 2. 
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Figure A9. The best x 2 for this source is z ~ 2.0, but secondary peaks are more consistent with its spectroscopic redshift (^spec — 
0.507). Longer wavelength data (e.g. at 850 /zm) are needed to constrain the redshift more accurately. The dust peak at ~ 200 fim. and 
the upper limits at longer wavelengths strongly constrain the source to be at z p hot < 3.8. 
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Figure All. A weak stellar bump peaking between the IRAC1 and IRAC2 bands is observed for this source. The overall SED is very 
similar to source 7 suggesting a low redshift. This is also found by the photometric redshift fitting procedure. The far-IR observations are 
not well fit but may be due to confusion with another source very bright at 24 jum and very close (~ 6 arcsec) to the center of detections. 
The IRAC photometry, however, is very constraining and the % 2 ditstribution also confirms a clear low redshift for this source. 
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